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Abstract 

The status of lattice calculations in Quantum Field Theory is reviewed. A major 
part is devoted to recent progress in formulating exact chiral symmetry on the 
lattice. Another topic which has received a lot of attention is the influence of 
dynamical quark effects. Attempts to quantify these effects for the light hadron 
spectrum and flavour singlet amplitudes arc discussed, as well as other expected 
qualitative features of simulations with dynamical quarks. The remaining parts of 
the review include recent results for the light quark masses using non-perturbative 
renormalisation and the spectrum of glueballs and heavy hybrids computed using 
anisotropic lattices. 
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1. Introduction 

Since Wilson's original formulation of lattice gauge 
theories in 1974 ||], lattice methods have developed 
into a mature area of research with a wide range of 
applications, including subjects as diverse as QCD, 
Higgs models and Quantum Gravity. It is therefore 
quite impossible to cover all recent activities in 
Lattice Gauge Theory in a single review. Therefore, 
I shall concentrate on lattice QCD applied to the 
calculation of hadron masses and the evaluation of 
weak matrix elements. Another subject which I 
shall present in detail, and which has received a 
lot of attention recently, is the lattice formulation 
of chiral symmetry. Indeed, over the past two years 
it has emerged that chiral gauge theories can be 
put on the lattice in a consistent way, something 
which has been thought to be impossible for a long 
time. These new developments have wide-ranging 
implications for a large class of theories, which 
clearly merits a detailed discussion. 

An overview of those areas which will not be 
discussed here can be found in the proceedings of 
recent annual conferences on Lattice Field Theory. 
In particular, I refer the reader to the plenary talks 
on 

• QCD at finite temperature and density (|-^ 

• The electroweak phase transition (5|-|^] 

• Topology and confinement [^,^ 

• Lattice gravity and random surfaces [|o[|ll| 

Furthermore there are reviews of lattice gauge 
theory given at recent conferences on High Energy 
Physics ||,|l|. 

t Plenary talk presented at International Europhysics Con- 
ference on High-Energy Physics (EPS-HEP 99), Tampere, 
Finland, 15-21 July 1999. 



1.1. General remarks 

Many of the concepts and techniques of Lattice 
Gauge Theories can be easily introduced in the 
context of QCD. Here, the lattice formulation 
provides a non-perturbative framework to compute 
relations between Standard Model parameters and 
experimental quantities from first principles. 

The discretisation is achieved by introducing 
a euclidean space-time lattice with spacing a and 
volume ■ T. The inverse lattice spacing, 
acts as an UV cutoff, which preserves the gauge 
invariance of the theory. The quark and antiquark 
fields ''Pix) are associated with the lattice 

sites X, whereas the gauge field is represented by 
the so-called link variable U^{x), which connects 
neighbouring lattice sites, and is an element of 
the gauge group SU(3). After choosing suitable 
discretisations of the Yang-Mills action Sq and 
the quark action 5f, the expectation value of an 
observable Vl is defined as 



1 



X,fJ, 



Y[dU^ix) Y[d^pix)di:{x)ne- 



Sq — Sf 



(1) 

and the functional integral Z is determined by 
requiring (1) = 1. The discretisation procedure has 
hence given a meaning to the functional integral 
measure, which becomes a simple product measure. 
More importantly, after integrating out the quarks 
this formulation allows for a stochastic evaluation 
of (il) using Monte Carlo techniques. 

Ideally one would like to use lattice QCD as 
a phenomenological tool. A typical application 
would be, for instance, the computation of the 
strange quark mass (2 GeV) , using the kaon 
mass, ruK and the pion decay constant as input 
parameters. However, as I shall explain below, 
realistic simulations of lattice QCD are difficult. 
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The first problem one has to address are lattice 
artefacts (cutoff effects). Let ft denote the quantity 
we wish to compute on the lattice, e.g. a hadron 
mass. Then the expectation values on the lattice 
and in the continuum differ by corrections of order 
aP, viz 

{nY''^ = {ny°''' + o{aP), (2) 

where the power p in the correction term depends 
on the chosen discretisation of the QCD action. 
Values of a which can currently be simulated lie 
in the range a « 0.2 — 0.05 fm. The size of the 
correction term can in some cases be as large as 
20%, depending on the quantity and the chosen 
discretisation. It is then clear that an extrapolation 
to the continuum, a — > is required in order to 
obtain the desired result. This extrapolation can be 
much better controlled if the chosen discretisation 
avoids small values of p. 

Perhaps the biggest challenge as far as lattice 
QCD is concerned, is the inclusion of dynamical 
quark effects. After integrating out the quark fields 
the expression for the expectation value reads 



-Sg 



where D is the lattice Dirac operator and my is 
the mass of quark flavour /. The evaluation of the 
determinant in eq. (|^) in numerical simulations is 
still very costly, even on today's massively parallel 
computers. In many applications the determinant 
has therefore been set to 1. This defines the so- 
called quenched approximation, which corresponds 
to neglecting quark loops in the evaluation of 
(fl). Although this represents a rather drastic 
assumption about the influence of quark-induced 
quantum effects, the quenched approximation works 
surprisingly well, as I shall describe later. 

An indirect consequence of using the quenched 
approximation is the observed scale ambiguity. 
That is, the calibration of the lattice spacing in 
physical units, [MeV], is dependent on the 
quantity Q which is used to set the scale 



'1 [MeV] 



Q [MeV] 
(aQ) ■ 



(4) 



This ambiguity arises because different quantities Q 
are affected by quark loops in different ways. 

There are also restrictions on the quark masses 
ruq that can be simulated. In general the following 
inequalities should be satisfied 



typical hadronic state, which serves as a measure 
of the quark mass. The inequality on the right 
places restrictions on the light quark masses that 
can be simulated: if those are too light one may 
suffer from finite-size effects, since ^ becomes large. 
Typical spatial extensions of i « 1.5 — 3fm imply 
that the physical pion mass cannot be reached. The 
left inequality restricts the masses of heavy quarks. 
Since a^^ w 2 — 4GeV, it is clear that relativistic b- 
quarks cannot be simulated. One therefore relies on 
extrapolations in ruq to connect to the physical w, d 
and b quarks. Obviously it is of great importance 
to control such extrapolations. 

Finally there is the issue of chiral symmetry 
breaking. A famous no-go theorem by Nielsen 
and Ninomiya [l^ ] implies that under fairly mild 
assumptions exact chiral symmetry cannot be 
realised at non-zero lattice spacing. Therefore, the 
chiral and continuum limits cannot be separated. 
For many (but not all) applications of lattice QCD 
this may not be a severe limitation, but the no-go 
theorem has so far precluded all attempts to achieve 
a realistic lattice formulation of the electroweak 
sector of the Standard Model. 

1.2. Outline 

The remainder of this article is as follows: 
section || deals with the properties of lattice Dirac 
operators: it is described how lattice artefacts 
in observables can be successfully reduced by 
constructing "improved" lattice actions. The 
major part of section ^ is devoted to recent 
developments which have ultimately led to the 
construction of lattice chiral gauge theories with 
exact gauge invariance. Section ^discusses recent 
simulations of QCD with dynamical quarks. Results 
which quantify sea quark effects in the light 
hadron spectrum are presented, and moreover the 
qualitative effects of sea quarks in flavour-singlet 
amplitudes and regarding the breakdown of linear 
confinement ("string breaking") are discussed. In 
section ^ recent results for the light quark masses 
are presented. Progress in this area has been 
achieved through the non-perturbative matching of 
lattice results to the MS scheme. In section E I 
describe recent results for the spectrum of glueballs 
and heavy hybrids. Section | contains a brief 
overview of results for some weak hadronic matrix 
elements which could not be reviewed extensively 
due to lack of time. Finally, a summary is presented 
in section 0. 



(|5^ 2. The lattice Dirac operator 



where L is the spatial extent of the lattice volume. 
The quantity ^ denotes the correlation length of a 



This section deals with the general properties of lat- 
tice fermions. After recalling the fermion doubling 
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problem we shall discuss the implementation of the 
Symanzik on-shell improvement programme, which 
systematically reduces lattice artefacts in physical 
observables and thus makes extrapolations to the 
continuum limit more reliable. The main part of 
this section, however, reviews the recent progress 
made in formulating chiral symmetry on the lattice. 

2.1. Fermion doubling revisited 

Suppose we want to describe massless free fermions 
on the lattice. The lattice action can be written as 

SY = a^Yj'{x)D{x-y)i:{y), (6) 

where D denotes the lattice Dirac operator. In 
particular, one would like to formulate the theory 
such that D satisfies the following conditions 

(a) D{x — y) is local 

(b) D{p) = ^7^p^ + Oiap") 

(c) D{p) is invertible for p 7^ 

(d) 75Z? + i?75 -0 

Locality is required in order to ensure renormalis- 
ability and universality of the continuum limit; it 
ensures that a consistent field theory is obtained. 
Furthermore, condition (c) ensures that no addi- 
tional poles occur at non-zero momentum. If this 
is not satisfied, as is the case for the "naive" dis- 
cretisation of the Dirac operator, additional poles 
corresponding to spurious fermion states can ap- 
pear: this is the famous fermion doubling problem. 
Finally, condition (d) implies that S-p is chirally in- 
variant. 

The main conclusion of the Nielsen-Ninomiya 
no-go theorem is that conditions (a)-(d) cannot 
be satisfied simultaneously. Since one is not 
willing to give up locality and condition (b), this 
implies that one is usually confronted with the 
choice of tolerating either doubler states or explicit 
chiral symmetry breaking. This is manifest in the 
two most widely used lattice fermion formulations: 
staggered ("Kogut-Susskind") fermions leave a 
chiral U(l) subgroup invariant, but only partially 
reduce the number of doubler species p5| . Wilson 
fermions, on the other hand, remove tEe doublers 
entirely at the expense of breaking chiral symmetry 
explicitly. This is easily seen from the expression 
for the free Wilson-Dirac operator 

<^ = H(VM + V;)-iaV*V^. (7) 

Using the definitions for the forward and backward 
lattice derivatives, and V*, one easily proves 
conditions (a)-(c), while it is obvious that (d) is not 



satisfied. Nevertheless, the (interacting) Wilson- 
Dirac operator Dw is widely used in simulations 
of lattice QCD. This is because chiral symmetry 
breaking in a vector-like theory is often merely an 
inconvenience, but no fundamental obstacle. 

2.2. Improved discretisations of the Wilson action 

Another consequence of using the Wilson-Dirac 
operator is the presence of large cutoff effects 
in physical observables. One can show that the 
leading lattice artefacts are of order a (i.e. p = 1 
in eq. (|^)), while e.g. for staggered fermions the 
lattice corrections start at order . 

It is important to realise the definition of 
a lattice action is not unique: one can add 
any number of operators which formally vanish 
in the continuum limit, provided that they 
comply with the correct symmetry and locality 
requirements. Given this relative freedom in 
defining lattice actions, an obvious question is 
whether one can find improved discretisations for 
lattice fermions, for which the cutoff effects are 
reduced. Two approaches have been the subject 
of much recent activity. One is based on the 
Symanzik improvement programme in 
which lattice artefacts are removed order by order 
in the lattice spacing. In the second approach 
one seeks to construct a "perfect" lattice action 
(which is essentially free of lattice artefacts) , using 
renormalisation group techniques . 

Sheikholeslami and Wohlcrt have shown that 
in QCD with Wilson fermions the Symanzik 
improvement programme can be implemented to 
lowest order in a by adding one dimension-5 
counterterm to Dw. The resulting operator (with 
bare mass mg) reads 

Dsw ^ Dw + ma + fcs^af,yFf,^, (8) 

where F^i, is a lattice transcription of the field 
tensor. In order to remove all lattice artefacts 
of order a in hadron masses, the improvement 
coefficient Csw has to be fixed by imposing a 
suitable improvement condition. Such a condition 
is provided by requiring that the restoration of 
the axial Ward identity holds up to terms of 
order Q. It has been appfied to determine Csw 
non-perturbatively for a large range of couplings for 
both quenched QCD |2| and also for rif = 2 
fiavours of dynamical quarks . 

In order to compute matrix elements of local 
composite operators such as vector and axial 
vector currents, one has to consider their correctly 
normalised and improved versions. For instance. 
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Figure 1. Scaling plot*^ for masses in the vector 
meson and nucleon channels for improved (circles) 
and unimproved (crosses) Wilson fermions [E2|. The 
linearity of the data computed with the improved action 
demonstrates the expected scaling behaviour. 



the renormalised axial current (A^)^ in the 0{a) 
improved theory reads 

(Ar)^ = Za(1 + bAaniq) {A^, + CAad^^P} (9) 

where 6a and ca are improvement coefficients, 
P is the pseudoscalar density, and Za is the 
renormahsation constant for the axial current. In 
refs. ||2l|,|2j] ca and Za have been determined non- 
perturbatively, and a general strategy to determine 
the improvement coefficients and renormahsation 
constants of all quark bilinears in the 0{a) improved 
theory has been presented in . Detailed studies 
of the matrix elements of such operators with much 
better controlled lattice artefacts should therefore 
soon be feasible. 

The remaining question is whether non- 
perturbative 0{a) improvement can be verified, in 
the sense that physical observables computed in the 
0{a) improved theory can be shown to approach 
their continuum limit with a rate proportional 
to Detailed scaling studies performed in 

the quenched approximation have shown that this 
is indeed the case for meson and baryon masses , 
as well as for matrix elements of vector and axial 
vector currents ||2^ . An example is shown in Fig. |^, 
taken from ref. P4 . 

In ref. [2^ the Symanzik improvement pro- 
gramme was extended to higher orders. The many 
relevant improvement coefficients have been deter- 
mined in (mean-field improved p9[) perturbation 



theory. 0(a) improvement has also been studied 
for the anisotropic Wilson action , for which the 
temporal and spatial lattice spacings are chosen dif- 
ferently. This offers an advantage in the computa- 
tion of heavy states, and we will return to this point 
in section ^. 

Another method to construct a lattice action 
with improved scaling behaviour was presented in 
ref. pyl. The idea is to use the operator Dsw 
in eq^^) on smoothed gauge configurations. Such 
configurations are obtained through a blocking 
procedure applied to spatial link variables as 
described in |32|. The resulting action is called 
the "fat-link" clover action. In a series of 
papers ||3^,0 it has been shown that fat-link clover 
actions exhibit good chiral properties: the additive 
renormahsation of the quark mass encountered for 
Wilson fermions is small (yet non- vanishing) , and 
the renormahsation factors Za and Zy of axial and 
vector currents are close to unity. Furthermore, 
scaling tests of various hadron masses have been 
carried out to test whether cutoff effects of order a 
have been eliminated. 

One concludes that improved discretisations 
of the Wilson actions are successful and play an 
ever more important role in calculations of the 
hadron spectrum and matrix elements. Their 
main advantage is that more accurate results 
in the continuum limit are obtained. It is 
important to keep in mind, though, that improved 
discretisations obtained through the Symanzik 
improvement programme and also the fat-link 
clover actions do not alleviate the problem of 
explicit chiral symmetry breaking. 

2.3. Exact chiral symmetry on the lattice 

A lot of progress has been made recently in the 
formulation of chiral fermions on the lattice. In fact, 
it has been shown how lattice chiral gauge theories 
can be formulated in a way which preserves locality 
and gauge invariance. Many new developments 
in this field have followed the rediscovery of the 
Ginsparg- Wilson relation |^5| 



(10) 



The significance of this relation had not been re- 
alised for a long time, since no non-trivial expres- 
sion for a Dirac operator D satisfying eq. ( p^ was 
known. It is remarkable that two constructions of 
such a solution have been developed independently 
at around the same time. 

The first goes back to Kaplan's proposal to re- 
alise chiral fermions in the domain wall fermion 
approach This was subsequently re-cast 

by Narayanan and Neuberger, who developed the 
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"overlap" representation of the chiral determi- 
nant |^,|3^]. Furthermore, Shamir and Furman ps| ] 
used a variant of Kaplan's work to derive the do- 
main wall fermion formulation for vector-like theo- 
ries like QCD The second solution to eq. ( |lC| ) 
was constructed in the perfect action approach men- 
tioned earlier. 

Here I will briefly discuss some of the concepts 
and applications. More details can be found in 
recent reviews on Domain Wall fermions pO| and 
the construction of lattice chiral gauge theories H, 
@- 

I shall first discuss the Domain Wall (DW) 
fermion formulation as described in ref. The 
basic idea, proposed in is to introduce an 
extra (fifth) dimension and to consider fermions 
coupled to a mass defect in the extra dimension. 
To make this more explicit, let x,y denote 4-dim. 
coordinates and s, s' the coordinates in the 5*^ 
dimension, which has finite length Ng. The gauge 
fields are trivial in the 5**^ direction, and the Dirac 
operator then has the form 

D'^^''{x,y)^D\x,y)5ss'+5{x-v)Di, (11) 

where D^^{x^y) is the usual Wilson-Dirac operator 
with a negative mass term, —M. The operator 
couples fermions in the extra dimension and 
contains the Dirac mass m. It can now be shown 
that for TO = and in the limit Ng oo there are 
no fermion doublers and, more importantly, chiral 
modes of opposite chirality are trapped in the 4- 
dim. domain walls at s = 1, N^. The physical, 
4-dim. fields defined by 

q{x) = P-Ripiix) + P-l'iI'nS^) 
q{x)^^r^^ix)P^+^,{x)PL (12) 

satisfy exact, continuum- like axial Ward identities 
at non-zero values of a for Ng —^ oo. This result, 
derived in formally establishes the correct 

chiral properties of the lattice regularised theory. 

In a real simulation one has to work at finite Ns 
so that the decoupling of chiral modes is not exact. 
It is possible to show, however, that the terms 
which break the chiral symmetry are exponentially 
suppressed. This point can be illustrated by 
discussing the relation between the pion mass and 
the quark mass in Domain Wall QCD. At finite a 
one has 

(aTO^)2 = C(aTO + ae-''^=) , (13) 

where C is a constant, m is the Dirac mass which 
appears in D-^, and 7 > 0. This expression shows 
that in the limit Ns —^ 00 the bare quark mass is 
only multiplicatively renormalised (i.e. there is no 
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Figure 2. The pion masd-jn the chiral limit plotted 
versus the size of the 5th dimension from ref. p3[ . 
Solid circles are the data obtained using the Wilson 
plaquette action for the gauge field, whereas open 
symbols denote data obtained using an improved gauge 
action (Iwasaki). The solid line is a fit assuming 
{am.Tr)'^ — A + B exp{—jNs), which yields a non- 
vanishing intercept of A = 0.048(9) in the limit Ns ^ 00 
(notation translates as: Ls —> Ns,mo M,mf — * m). 

additive quark mass renormalisation as for ordinary 
Wilson fermions). It also shows that for finite a 
there is a residual pion mass in the chiral limit 
proportional to aexp(— 7iVs). 

Hence, the Domain Wall formulation of QCD 
offers a method to realise almost exact chiral 
invariance at non-zero lattice spacing at the expense 
of simulating a 5-dim. theory. An important feature 
is that lattice artefacts of order a are exponentially 
suppressed at finite Ng. In other words, the theory 
is 0{a) improved for iVs — s- 00, with exponentially 
small 0{a) corrections expected at finite Ns- As 
in the case of non-perturbative 0{a) improvement 
the expected improved scaling behaviour in Domain 
Wall QCD must be verified. 

The question that arises is how small a value 
of Ns one can get away with in order to realise 
chiral symmetry whilst keeping the computational 
overhead of simulating an extra dimension at a 
minimum. In refs. |E3,H4] systematic studies of the 
A^s-dependence of (otoT)'^ in the chiral limit have 
been presented. An example of such an analysis is 
shown in Fig. || |^^. The main conclusion is that 
the suppression of the residual pion mass expected 
according to eq. ( p^ is very slow. At present it 
is not even clear that it vanishes at all in the 
limit As — > 00. Further studies, clarifying the role 
of conventional, 3-dim. finite-size effects and the 
effects of quenching, are required to settle this issue. 

There are already extensive numerical studies 
of physical observables using DWQCD. Results 
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include studies of QCD thermodynamics , weak 
interaction matrix elements such as those relevant 
for kaon physics, i.e. the _B-parameter Bk and 
e'/e p6|-^ an d determinations of the light quark 
masses [H9fEo|. 

Turning now to the question of constructing 
chiral gauge theories on the lattice with exact gauge 
invariance, let us recall that a central ingredient in 
the construction of such theories is the Ginsparg- 
Wilson (GW) relation eq. (Oh , which replaces the 
condition _{75, D} = 0. (Jne key observation, 
made in |51t| , is that any lattice Dirac operator 
satisfying the GW relation has chiral zero modes 
and satisfies an exact index theorem. It was 
subsequently realised |Q that the GW relation 
implies an exact symmetry of the associated action, 
with infinitesimal variations proportional to 



J^ = 75(l-iaD)i/, 



(14) 



Moreover, this symmetry reproduces the correct 
chiral anomaly in the flavour singlet case. That is, 
all the hallmarks of the correct chiral behaviour are 
present in the lattice theory: chiral zero modes, an 
exact index theorem and the chiral anomaly derived 
from the Ward identities associated with the exact 
symmetry. 

At this stage left- and right-handed fermions 
are easily introduced. Furthermore, starting from 
a solution to the GW relation, it is possible 
to construct abelian chiral gauge theories on the 
lattice, which comply with all basic requirements, 
including exact gauge invariance [ |53| . The 
construction extends to the non-abelian case 
and although a few properties have yet to be 
established with the same rigour as for abelian 
theories, there is little doubt that the construction 
is valid. 

An operator which satisfies the GW relation 
can be constructed in the framework of the 
overlap formalism. After Kaplan's proposal 
and independent related work by Frolov and 
Slavnov pq |, Narayanan and Neuberger obtained 
an expression for the chiral determinant, which is 
usually termed the "overlap" [p7 , 38 1 . Among many 
other things the overlap was shown to reproduce 
the correct chiral properties in vector-like theories 
like QCD Furthermore, it can be brought 

into a more manageable expression in terms of a 
relatively simple lattice Dirac operator D^. It is 
defined by pi] 



D^ = Ul-X{X^X) 



-1/2 



where Dw is the massless Wilson-Dirac operator. 
_Dn can easily be shown to satisfy the GW 
relation ||5^. Moreover, it was demonstrated that 
the construction of Dn from the overlap could 
be reversed, using only the GW relation and 
factorisation Q. 

The obvious question at this point is whether 
the correct chiral properties of operators satisfying 
the GW relation can be verified in lattice 
simulations. A possible strategy is then to 
check whether the expected global anomalies are 
exhibited. One example is Witten's observation 
that SU(2) gauge theory coupled to a single left- 
handed fermion is mathematically inconsistent. 
Restricting the discussion to the continuum theory 
for the moment, this is easily seen by examining 
the functional integral for this theory after the Weyl 
fermions have been integrated out, viz 



Z = y [dA] {dct D[A]) 



1/2 q-Sg[A]^ 



(16) 



Here, A is the SU(2) gauge field, and D denotes the 
Weyl operator. Obviously a sign ambiguity arises 
due to the presence of the square root in eq. (^6|). 
Indeed, for SU(2) there exist non-trivial gauge 
transformations g{x) which cannot be deformed 
continuously to the identity, and for which 

(det D[A])^/^ = - (det D[A3])^^^ , (17) 

where A^ denotes the transformed gauge field. 
This result implies that expectation values defined 
through Z are indeterminate, i.e. (O) ="0/0". 

Although A^ and A^ are not connected via some 
smooth gauge transformation, they are none the 
less connected in the space of gauge fields. Thus, 
one can define a curve which smoothly interpolates 
between and AJ^, given by 



A^it) = {l-t)A^ + tAl te[0,l]. 



(18) 



X 



D 



(15) 



By invoking the Atiyah-Singer index theorem, 
Witten showed that the number of eigenvalues of 
the Dirac operator which cross zero is odd. This 
observation about the behaviour of the spectral 
flow along the curve A^{t) then gives rise to 
the sign ambiguity and hence the mathematical 
inconsistency. It is now interesting to investigate 
whether the expected spectral flow can be 
reproduced on the lattice. 

Initially this has been investigated by Neu- 
berger |^0|. Then, in ref. Bar and Campos re- 
ported on the computation of the eigenvalues of the 
overlap operator Dn in a lattice simulation along 
the path which smoothly connects a constant gauge 
configuration to its gauge transform as in eq. (ISh. 
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Figure 3. Level crossing for the imaginary part of the 
lowest eigenvalue of the overlap operator 



Focussing on the six lowest eigenvalues they found 
that only one, i.e. the smallest, Aq, becomes zero 
at t = 0.5. The hermiticity properties of im- 
ply that one expects the level crossing to occur for 
the imaginary part of this eigenvalue. Figure 
demonstrates that the crossing is indeed observed 
for ImAo. This result represents a numerical proof 
of Witten's original argument and is only possible 
if the lattice Dirac operator has the correct chiral 
properties. The above example illustrates the enor- 
mous progress that has been achieved in the formu- 
lation of chiral symmetry. 

Other recent work in this area includes detailed 
studies of the axial anomaly ]6^-|67[|, investigations 
of spontaneous chiral symmetry breaking [6^ - 
[TOI , the locality properties of Z?n [0 and 
the development of efficient implementations of 
Ginsparg- Wilson fermions ||7^-[7^ . 

In view of these results it should be clear 
now that a consistent formulation of the Standard 
Model exists beyond perturbation theory. There 
are also important consequences for future lattice 
studies of supersymmetric models: in many ways 
the intrinsically supersymmetric features in such 
models bear resemblance to the role of chiral 
symmetry in QCD. 

3. Simulations with dynamical quarks 

Now we turn to discuss the other great challenge in 
present simulations of QCD, namely the inclusion 
of dynami cal quark effects. As mentioned 
in section 
still widely 



1.1 



the quenched approximation is 
used for many phenomenologically 
interesting quantities, and in order to enhance the 
predictive power of lattice simulations it is of great 
importance to assess the influence of dynamical 
quarks. 



Before such an assessment can be made with the 
required accuracy, it is necessary to perform precise 
calculations of experimentally known quantities in 
the quenched theory. Such benchmarks serve not 
only to detect significant effects due to dynamical 
quarks, but also illustrate how closely the quenched 
approximation resembles the real world. 

3.1. Quenched light hadron spectrum 

The CP-PACS Collaboration recently presented a 
precision calculation of the light hadron spectrum 
in quenched QCD using the (unimproved) Wilson 
action 
GFll 



which superseded earlier studies by 
The findings of CP-PACS are 
summarised in the plot shown in Fig. |[ Although 
the qualitative features of the spectrum are well 
reproduced by the quenched lattice data, one 
finds significant deviations from the experimentally 
observed spectrum. For instance, the ratio of the 
nucleon and p masses is calculated as 



mN/mp = 1.143 ± 0.033, 



(19) 



which is 6.7% (2.5 fi) below the experimental 
value of 1.218. Similarly, vector-pseudoscalar mass 
splittings such as ttik* — itt-k are too small by 
10—16% (4 — 6(t), depending on whether mx 
or is used as input to fix the strange quark 
mass. This result implies that, for the first 
time, a significant deviation between the quenched 
QCD spectrum and nature is detected. The 
main conclusion of ref. is that quenched QCD 
describes the light hadron spectrum at the level 
of 10%. However, it also shows that the quenched 
approximation works surprisingly well, since the 
discrepancy is fairly mild. This has important 
consequences for other, phenomenologically more 
interesting quantities, for which one may be stuck 
with the quenched approximation for some time to 
come. 

Although the CP-PACS results represent a real 
benchmark in terms of statistics, parameter values 
and lattice volumes, further corroboration of these 
findings is still required. It should be added that 
the results depend crucially on the modelling of 
the quark mass dependence of observables. Usually 
the results of Chiral Perturbation Theory can be 
used to guide the extrapolations in the quark 
masses. However, in the quenched approximation 
one expects deviations from the predictions of 
Chiral Perturbation Theory due to the appearance 
of "quenched chiral logarithms" ^l|, which 
makes the extrapolation of hadron masses close to 
the chiral limit hard to control. Furthermore, the 
CP-PACS results were obtained using unimproved 
Wilson fermions which have large lattice artefacts. 
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Figure 4. The quenched hadron spectrum from ref. j?^ 
compared to experiment (dashed hues). 



It is therefore desirable to check whether the 
extrapolations to both the chiral and continuum 
limits are controlled. 

Recent calculations employing different discreti- 
sations have largely confirmed the findings of [jTZf : 
the MILC Collaboration |82| has used staggered 
fermions and finds a value for the nucleon to rho 
mass ratio in the continuum limit of 

mjv/mp = 1.254 ± 0.018 (stat) ± 0.028 (syst), (20) 

This is in broad agreement with experiment, but 
the difference to the CP-PACS result amounts to 
only two a. A recent calculations by UKQCD using 
0(a) improved Wilson fermions |83] essentially 
confirms the conclusion of CP-PACS7namely that 
the quenched light hadron spectrum agrees with 
experiment at the level of 10%. 

3.2. Light hadron spectrum for n[ = 2 

An obvious question is whether sea quark effects can 
account for the observed deviation of the quenched 
light hadron spectrum from experiment. Before 
discussing some results it is useful to explain a few 
technical issues which are relevant for simulations 
with dynamical quarks. 

First it is evident that improved lattice 
actions have an even more important role to 
play. Since simulations with dynamical quarks 
are very expensive, one cannot afford to control 
the extrapolation to the continuum limit by 
simulating very small lattice spacings, whilst 
keeping sufficiently large spatial volumes in physical 
units. In order to be able to separate sea 
quark effects from lattice artefacts, it is vital to 



Table 1. Recent simulations with rif = 2 flavours 
of dynamical quarks, together with the sustained 
computer power (in GFlops). The other two columns 
denote the choice of discretisation for the gauge action 
So and the Dirac operator used in Sf ■ 



CoUab. 


GFlops 


Sg 


S'f 


RBC 


- 250 


Plaq. 


£)DWF 


CP-PACS 


- 300 


Iwasaki 


Dsw, tad. 


UKQCD 


- 30 


Plaq. 


Dsw, n.p. 


SESAM/TxL 


- 14 


Plaq. 


Dw 


MILC 


> 7 


LW 


^stag 



have control over the latter. Details of recent 
simulations with n{ = 2 flavours of dynamical 
quarks are listed in Table |l]. The quark actions 
used in these calculations are based on the 0{a) 
improved Wilson-Dirac operator Dsw (using either 
the non-perturbative determination for Cgw or 
its estimate in mean-field improved perturbation 
theory) , the unimproved operator Dw , domain wall 
quarks or staggered fermions. In addition to the 
widely used plaquette action for the gauge fields 
(labelled "Plaq." in Table |l|), two collaborations 
also use improved gauge actions. Such actions 
were introduced by Liischer and Weisz |Q and 
Iwasaki they have been found to lead 

to smaller lattice artefacts in the static quark 
potential, although the bulk of the discretisation 
effects in hadronic quantities can only be reduced 
by employing improved fermion actions |]86|| . 

The two flavours of (degenerate) dynamical 
quarks with mass m^°'^ are usually identified with 
the physical u and d quarks. Unlike in the 
real world it is possible in lattice simulations to 
compute observables relating to hadrons whose 
valence quarks have a different mass than the sea 
quarks, i.e. to™' ^ m^°^. This is illustrated in the 
diagram in Fig |^ for the case of a mesonic two-point 
function. It implies that one has more freedom 
to explore the dependence of physical observables 
on TO™' and m*"^^ separately, and to compare the 
results to the predictions to Chiral Perturbation 
Theory ||79| and its modified version for unequal sea 
and valence quark masses |^-^ . 

It is clear that the simulations listed in Table |^ 
do not have the correct value of rif for kaon physics. 
In the absence of any sufficiently tested simulation 
algorithm capable of treating rif = 3 flavours of 
dynamical fermions one is forced to introduce the 
strange quark as a valence quark by making the 
identifications 



^ val 



mu,d 



(21) 
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Figure 5. Quark and gluon contributions to the two- 
point function for a pseudoscalar (vector) meson in 
partially quenched QCD with m^^" / rn""^. 



This defines the so-called partially quenched 
approximation. 

One should also bear in mind that the values 
of m^°^ which are accessible in current simulations 
(especially for Wilson quarks) are still relatively 
large. This is easily seen by computing the ratio 
of the pseudoscalar to vector meson mass mps / my 
(for rrf^^ = m™^) and comparing it to the physical 
ratio mT^/rnp = 0.169. For some of the simulations 
in Table | one finds @,|9^-|9| 



my 
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0.58 
0.60 



0.83 
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0.80 
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CP-PACS 



(22) 



In order to make contact with the physical situation 
one has to study the dependence of observables on 
^soa extrapolate in m^'^^ to the physical value of 
rriT^/mp. Since the vector meson can decay into two 
pseudoscalars below mps/rriY « 0.5 such a naive 
extrapolation may, however, be misleading. 

Several collaborations have investigated sea 
quark effects in the light hadron spectrum. CP- 
PACS studied the continuum limit of hadron 
masses computed for rif = 2 and compared it with 
the results of the quenched light hadron spectrum 
discussed before. Figure ^ shows the results for the 
masses of the K* and mesons for nf = 0, 2. This 
demonstrates clearly that the discrepancy with the 
experimentally observed spectrum is considerably 
reduced when sea quarks are "switched on" . The 
figure also illustrates that a reliable quantification 
of sea quark effects can only be made after the 
extrapolation to the continuum limit: at non-zero 
values of a the deviation between experiment and 
partially quenched QCD is enhanced, and thus one 
would come to the wrong conclusion about the size 
of dynamical quark effects. 

Results for the vector-pseudoscalar mass split- 
ting reported by UKQCD |91| also show that lattice 
data for this quantity approach the experimental 
value as the sea quark mass is decreased. 

These findings are quite encouraging: they 
demonstrate that sea quarks have the expected 
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Figure 6. Scaling behaviour of vector mesons in 
quenched and partially quenched QCD [ p3[ . 



effects on hadronic quantities. Of course, the 
remaining differences between partially quenched 
QCD and experiment seen in Fig. ^ have to be 
explained. It is reasonable to assume that the 
incorrect value of nf for kaon physics will have some 
influence. Further studies are required in order to 
decide whether there is yet sufficient control over 
the extrapolations in m^'^'^, to™', as well as the 
extrapolations to the continuum limit. 



3.3. Sea quark effects in other observables 

There are a number of observables which are 
particularly sensitive to sea quark effects. One 
example are flavour-singlet amplitudes which are 
relevant for quantities such as the 77' mass. 
Unlike flavour non-singlet matrix elements these 
amplitudes receive contributions from disconnected 
diagrams and are thus sensitive to vacuum 
polarisation effects due to dynamical quarks and the 
value of Uf. 

In a number of recent papers calculations 
of disconnected diagrams contributing to the rj' 
mass |9^ , the 7r-nucleon cr-term |9^ and the 
flavour-singlet axial coupling of the proton j9^ have 
been presented. Although some of the expected 
qualitative features of sea quark contributions 
to these observables have been seen the 
main obstacles for more quantitative analyses are 
the high level of statistical noise encountered 
in the evaluation of disconnected diagrams and 
also the fact that nf = 3 flavours cannot be 
simulated. It is estimated that the calculation 
of the disconnected contribution to the mass of 
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Figure 7. The static qu£iJc '"potential and the first 
excitation computed in th e SU(2) Higgs model using a 
variational technique [ |l06[ . The flattening of the ground 
state for large separations indicates the breakdown of 
linear confinement. 



flavour-singlet mesons is an order of magnitude 
more expensive than conventional correlators, so 
that more efficient numerical methods have to be 
developed and applied |98|. Therefore, there are no 
firm quantitative results for vacuum polarisation in 
flavour-singlet amplitudes at present. The current 
status of this field has been reviewed extensively 
in M]. 

A phenomenon which is expected to occur in 
the presence of dynamical quarks is the breakdown 
of linear confinement, also called "string breaking" . 
This should manifest itself in a fiattening-off of the 
static quark potential at a characteristic distance 
at which the mass of two heavy-light mesons is 
energetically favoured over the energy of the fiux 
tube which is responsible for the linearly rising 
potential between the static quarks. In QCD 
string breaking can be lin ked to decay rates of 
processes like T(4S) BB [100| and is therefore of 
direct phenomenological relevance. Despite many 



recent efforts [^6 
sign for string 



|91|,|101 there is no unambiguous 
Drcaking in simulations of QCD 
with rif = 2 fiavours of Wilson fermions, if the 
potential is determined by measuring Wilson loops. 
It has, however, bee r i ob served for QCD with 
staggered qu arks |102, 103(1 and in QCD at finite 
temperature [I04|. 

The reason for the failure to detect string 
breaking in QCD at zero temperature has been 
attributed to the bad projection properties of 
Wilson loops onto the state of the broken string. 



This has been confirmed in studies in which the 
fermionic fields have been replaced by scalar fields 
(which are computationally much less demanding 
whilst preserving the under lying m echanism for 
string breaking to occur) |105-107|. In these 
studies string breaking was treated as a mixing 
phenomenon between the string and the two- 
"meson" state. This was achieved by supplementing 
the operator basis with operators having an explicit 
projection onto the broken string and by employing 
a variational approach to determine the energy 
levels. Indeed this more sophisticated method has 
provided clear evidence for string breaking (an 
exarnple for the SU(2) Higgs model is s^hown in 
Fig. Zf), and there are preliminary results |108 , 109 1 
which indicate that the approach is successful in 
QCD as well. 

4. Light quark masses 

Quark masses are fundamental parameters of the 
Standard Model, and their determination has been 
the subject of many recent activities. In particular, 
the mass of the strange quark, nzs, has received a 
lot of attention foll owing the recent experimental 
results on e'/e | |110( |. Theoretical analyses of this 
quantity rely on mg as one of the essential input 
parameters. 

Ratios of the light quark masses are predicted 
by Chiral Pert urba tion Theory with a precision of 
a few percent [111|. In order to obtain individual 
quark masses it is thus sufficient to determine a 
particular linear combination using lattice QCD. 

^.1. Non-perturbative renormalisation 

A convenient starting point to discuss lattice 
calculations of light quark masses is the PC AC 
relation. For charged kaons it can be written as 



/k"1k = ("T-u + ms)(0|u75s|K). 



(23) 



In order to determine the sum of quark masses 
(tou -\-ms) using the experimental result for /kti^, 
it suffices to compute the matrix element (0|u75s|K) 
in a lattice simulation. The renormalisation of the 
pseudoscalar density M75S is, however, scale- and 
scheme-dependent. By convention quark masses 
are quoted in the MS scheme of dimensional 
regularisation at some reference scale ji. Therefore 
the relation between the pseudoscalar densities in 
the MS scheme and lattice regularisation must be 
computed, viz 



(m75s)ms = Zp{gQ,afi) {wy^s) 



lat 



(24) 

Here go is the bare coupling, and /i is the subtrac- 
tion point in the MS scheme. The renormalisation 
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factor Zp has been computed in lattice perturba- 
tion theory for several discretisations. However, the 
limitations of lattice perturbation theory are well 
known, and in order to remove all doubts about the 
reliability of the matching procedure it is evident 
that a non-perturbative determination of the renor- 
malisation factor is required. 

One of the main obstacles for a fully non- 
perturbative matching procedure are the large scale 
differences between the low-energy regime, which 
is the domain of the lattice regularisation scheme 
and the high-energy, perturbative regime where the 
MS scheme is defined. This technical difficulty 
can be overcome by introducing an intermediate 
scheme X, as shown schematically in Fig. ^. The 
problem is thus split into two parts. The first is the 
matching of the bare current quark mass m.iat(a) 
to the running mass in the intermediate scheme, 
rnxip)- This amounts to computing Zp{go, afio) 
between the lattice scheme and scheme X for a 
range of bare couplings go at a fixed scale /ip. 
The second part is the determination of the scale 
dependence of the running mass mxip-) from fio 
up to very high energies, where the perturbative 
relation between mx and fn-y^ can be expected to 
be reliable. Through this two-step process the use 
of lattice perturbation theory is completely avoided. 

So far two proposals to define a suitable 
intermediate scheme have been put forward. The 
first has been introduced in ref. [112| and imposes 
non-perturbative renormalisation conditions on 
Green functions of local operators, computed 
between off-shell quark and gluon states, with 
virtualities /i, in a fixed gauge. In the case of the 
pseudoscalar density such a normalisation condition 
reads 



Zp{go, a^i) Z ^{go, a^i) Tp{ap) 



1. 



(25) 



Here Fp denotes the amputated Green function 
of 975(7, and Zq is the quark wavefunction 
renormalisation factor (which is easily computed 
in a lattice simulation). Typically this relation 
is evaluated in the Landau gauge. This choice 
of intermediate scheme is referred to as the 
Regularisation Independent (RI) scheme, in which 
the role of the renormalisation scale is played 
by the virtuality of the quark states. The 
scale dependence can then be probed by choosing 
different external momenta for the quark fields used 
to compute Fp. Provided that the virtualities can 
be fixed such that 



Aqcd < M < a 



(26) 



the perturbative matching between the RI and MS 
schemes expected to work well. 
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Figure 9. Non-perturbative scale evolution of msp/M 
computed in lattice simulations of the SF (solid circles). 
The dotted, dashed and solid lines correspond to the 
scale evolution computed using the 2/1-, 2/2 and 
3/2-loop expressions for the KG /3- function and the 
anomalous dimension. 



Another intermediate scheme is defined using 
the Schrodinger Functional (SF) of QCD p^|lT5|| . 
This scheme is based on the formulation of QCD in 
a finite volume of size ■ T with inhomogeneous 
(Dirichlet) boundary conditions in the time direc- 
tion. Non-perturbative renormalisation conditions 
can then be imposed at scale /i = 1/L and zero 
quark mass. An attractive feature of the SF scheme 
is that is allows to compute the scale dependence 
non-perturbatively over several orders of magni- 
tude, using a recursive finite-size scaling technique. 
Thus, once the scale dependence of msF is known up 
to energies of around 100 GeV one can continue the 
scale evolution to infinite energy using the pertur- 
bative renormalisation group functions and thereby 
extract the renormalisation group invariant (RGI) 
quark mass M. In Fig. ^ the scale depe ndence of 
TOgp/M computed in quenched QCD |116] is shown 
and compared to the perturbative scale evolution. 

The left-most point in Fig. ^ corres pond s to a 
scale ^0 ~ 275 MeV. Here one reads off |116[ 



M 
msF 



= 1.157 ±0.015. 



(27) 



The matching between lattice regularisation and 
MS scheme via the SF is completed by computing 
the renormalisation factor Zp{go^afj,o) for a 

particular fermion discretisation, at fixed /iq = 

275 MeV, for a range of bare couplings. In ref. [116| 
this has been performed for the 0(a) improved 
Wilson action, for couplings which correspond to 
lattice spacings in the range 



0.1 - 0.045 fm. 



(28) 
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Figure 8. Schematic relation between quark masses in lattice regularisation and the MS scheme through an 
intermediate renormalisation scheme X. 



4-.2. Light quark masses in quenched QCD 

Estimates for the light quark masses for several 
fermion lattice actions and using non-perturbative 
renormalisation (implemented either in th e RI or SF 
schemes) have been reported recently ^0, 117 119 1. 



As an illustration I shall describe in some detail 
the determinati on o f the strange quark mass along 
the lines of ref . |118[ , in which the non-perturbative 
renormalisation factors computed in ref. 116 were 
used. 

The average light quark mass mi is defined as 



mi 



i(mu 



md) 



By combining the result for M/mgp from eq. ( |27| ) 
with the factor Zp ((70,0/^0) a-nd the PCAC 



relation eq. (23) one obtains an expression for the 
sum of RGI quark masses (Mg -I- Mi) in units of the 
kaon decay constant 



Ms + Ml _ M ml^ 

/k msF Zp{go,afj.o)GKia) 



+0{a^). (30) 



Here, up to a small mass dependent factor which 
arises in the 0{a) improved Wilson theory, Gk is 
equal to the matrix element of the pseudoscalar 
density evaluated at the kaon mass. At this point all 
reference to the intermediate SF scheme is cancelled 
in the product M/(tosf^p)- What remains to be 
specified is the experimental value of expressed 
in units of some quantity which sets the lattice 
scale (e.g. the hadronic radius rg |121, 122|). The 
results for (A/g + Mi)/ as obtained from eq. ( |3C| ) 
can now be extrapolated to the continuum limit. 
This is shown in Fig. |l^, which also illustrates 
that the remaining discretisation errors in the 
0(a) improved theory are consistent with a leading 
behaviour proportional to a^. 

Using the experimental value /k = 160 ± 2 MeV 
one obtains in the continuum limit 



Ms + Ml = 140 ± 5 MeV. 



(31) 



This result can now be converted into m^®(/i) at 
/i = 2 GeV. First one combines eq. (|3^) wi th th e 
prediction from Chiral Perturbation Theory [ [lll[ 

Ms/Ml = 24.4 ± 1.5. 
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Figure 10. Continuum extrapolation of (Ms + Afi)//K. 
The lattice scale is set by hadronic radius ro. 



By integrating the (4-loop) perturbative RG 
functions in the MS scheme one obtains 

m'^(^)/M = 0.7208 at n^2GeY. (33) 

Finally, combining eqs. (^, ( ^ and (|3^) yie lds t he 
final result in the quenched approximation ]118| 



(2 GeV) = 97 ± 4 MeV. 



(34) 



(32) 



The quoted uncertainty of ±4 MeV contains all 
errors, except those due to quenching. As 
mentioned in the introduction, the conversion 
into physical units is ambiguous in the quenched 

approximation. For r?T^^(2GeV) the resulting 
uncertainty was estimated to amount to 10%. 

A compilation of recent results for TOs and mi 
in quenched QCD is shown in Table |^. When 
comparing the results one has to bear in mind 
that systematic errors - where shown - have not 
been estimated in a uniform manner, or have 
sometimes been combined in quadrature with the 
statistical errors. Also, the conversion into physical 
units has been performed using different quantities. 
Nevertheless, the picture that emerges is quite 
encouraging. Estimates for the strange quark mass 
in quenched QCD cluster around 100 MeV and 
around 4.5 MeV for the average up and down quark 
mass. Different discretisations such as Wilson, 
staggered and Domain Wall fermions as well 
as different implementations of non-perturbative 
quark mass renormalisation yield broadly consistent 
results. 
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Table 2. Estimates for the quark masses ttIs and mi in the MS scheme a,t ^ = 2 GeV in quenched QCD. Also 
shown is the choice of intermediate renormalisation scheme X and the lattice Dirac operator. Crosses indicate 
that non-perturbative renormalisation and/or a continuum extrapolation has not been implemented. 
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Figure 11. The strange quark mass and the ratio 
ma /mi extracted using either the vector or axial vector 
Ward identities plotted versus the lattice spacing. 



4-3. Sea quark effects in mg and m\ 

An important issue, especially for phenomenological 
applications of lattice calculations of the light 
quark masses is the influence of dynamical quarks. 
Estimates for mi and rrig computed in p art ially 



quenched QCD have been reported in |90, 93, 125 
|12q]. Since non-perturbative renormalisation has 
so far not been applied for rif = 2, all of these 
calculations rely on perturbation theory to match 
the quark masses in the lattice and MS schemes. 

The most comprehensive study so far has been 
presented by CP-PACS in which not only the 
m'^'^'^-dependence has been investigated but also 



the extrapolation to the continuum limit. The 

scaling behaviour of fn^^ (2 GeV) extracted either 
from the vector or axial vector Ward identities 
(labelled VWI or AWI respectively) is shown in 



the upper part of Fig. 11. In spite of the large 



differences observed for the two methods at non- 
zero lattice spacing, the results are consistent with 
a common continuum limit, where TO^^(2GeV) = 
84 ± 7MeV. This is significantly lower than 
the quenched results discussed earlier. However, 
many systematic effects are not as well controlled 
as in the quenched case, so that further studies, 
employing non-perturbative renormalisation and 
investigating the rif-dependence are required before 
a substantial decrease of quark masses relative 
to the quenched results can be confirmed. It is 
interesting, though, that the ratio mg/mi, in which 
many systematic effects are expected to cancel, 
extrapolates to ms/m\ = 26 ± 3 (c.f. lower part 
of Fig. ^l]). This value is in excellent agreement 
with the prediction from Chiral Perturbation 
Theory, eq. (|3^). T he pr eliminary results by the 
MILC Collaboration |126], which are also shown in 
the figure, have been computed at a fixed value 
of m^'^'^ corresponding to m.„/mp = 0.56. This 
may explain why they are different from the CP- 
PACS results, which in turn have been extrapolated 
in nf^^ to the physical value. 

For phenomenological applications one may be 
tempted to convert the results presented here into 
a global estimate. Before more thorough studies of 
dynamical quark effects become available, I consider 
the quenched results (e.g. eq. (0)) the most reliable 
estimate. By accounting for the systematic errors 
due to using the quenched approximation one can 
then quote a global result for mg as 



(2 GeV) = 100 ± 5 (stat) (syst) MeV. (35) 



14 



The systematic error in the above estimate 
incorporates the afore-mentioned scale ambiguity 
of lOMeV, as well as the observed decrease in the 
central value for Uf — 2 dynamical flavours. 

5. Glueballs and heavy hybrids 

Historically glueball masses were among the first 
quantities to be computed in lattice gauge theories. 
Over the years the calculations have become more 
refined, for instance, by constructing efficient 
glueball operators, by studying the approach to 
the continuum limit, and by including higher spin 
states. 

The main difficulty in glueball calculations is 
the relatively high level of statistical noise in the 
correlation functions from which the masses are 
extracted. Here I shall focus on recent calculations 
in which this problem has been alleviated by 
using anisotropic lattice actions. Hence, from 
now on I shall distinguish between spatial and 
temporal lattice spacings, as and at, respectively. If 
Cg (t) denotes the correlation function of a glueball 
operator G{x), then its asymptotic behaviour for 
large separations t is given by 



CGW = 5](G(x,t)Gt(0)) 



-(atMG){t/at) 



(36) 



That is, the exponential decay of Ccit) is governed 
by the glueball mass in units of the temporal lattice 
spacing. Therefore, the larger atMQ, the quicker 
the decay of Cg (t) , so that its asymptotic behaviour 
may be difficult to isolate before the statistical noise 
becomes too large. By introducing an anisotropic 
lattice action with at <C , one can simultaneously 
achieve slow exponential fall-off of Ccit) whilst 
preserving large spatial volumes in physical units. 
The spatial volume in lattice units, however, can be 
kept small if Us is large in physical units, so that the 
calculations are more manageable. Typical values 
of as are 

as w0.2-0.4fm, (37) 

while the "aspect ratio" ^ = as /at is usually taken 
in the range 

^ = 3-5. (38) 

The idea of using anisotropic lattices in glueball 
calculation s is not new: it was used already 16 
years ago |127], before the advent of "smearing" 
techniques to construct glueball operators for 
which the asymptotic behaviour in the correlation 
function sets in very quickly. 

The idea of anisotropic lattice actions has been 
revived by the desire to compute the higher glueball 
states more accurately and to use very coarse spatial 



lattices so that these calculations can be performed 
on smaller computers. However, by pushing 
to very large values of a^ one may suffer from 
uncontrollably large discretisation effects. Here, 
the Symanzik improvement programme can be used 
to eliminate the leading lattice artefacts, so that 
the extrapolation of results obtained on coarse, 
anisotropic lattices is sufficiently controlled. 

In the rest of this section I shall mainly discuss 
recent results for the spectra of glueballs and heavy 
hybrids obtained using anisotropic actions. For 
more compreh ensive reviews I refer the reader to 
refs. |l2|,|29|. 



5.1. Glueball spectrum in quenched QCD 



In two recent papers Morningstar and Peardon [130| 
have presented results for the quenched glue- 
ball spectrum below 4 GeV computed using an 
anisotropic, 0{a^) improved gluon lattice action 
with aspect ratios ^ = 3 and 5. Thus, improve- 
ment has been employed to reduce artefacts asso- 
ciated with the spatial lattice spacing only. Non- 
perturbative determinations of the relevant im- 
provement coefficients are not available for the case 
at hand, so that one relies on their estimates in 
mean-field improved perturbation theory. This pro- 
cedure does not completely remove lattice artefacts 



of order aj. 
of order 



One thus expects leading cutoff effects 



2 2 

9 as, 



al 



(39) 

Since the lattice breaks rotational symmetry, 
glueball operators are, as usual, constructed 
from representations of the octahedral group: 
Ai,A2,E,Ti,T2. Figure |l| shows the scaling 
behaviour of the glueball masses in the PC = ++ 
channel, extracted from the various representations. 
The plot nicely demonstrates the restoration of 
rotational symmetry in the limit a^ — s- 0. For 
instance, both the E and T2 representations 
describe the spin-2 glueball in the continuum limit, 
but differ significantly at large, non-zero lattice 
spacings. 

The observed curvature in the continuum 
extrapolation of the scalar (0++ and 0*"*^'*') glueballs 
suggests that both 0{al) and 0(af ) lattice artefacts 
are sizeable for these states. Hence, if lattice 
spacings as large as as « 0.4 fm are included one has 
to use a more complicated model function for the 
continuum extrapolation, which ultimately leads to 
a loss of statistical precision for these states. 

The fin al, c ontinuum results for glueball states 
from ref. |130| for a number of different J^'-^ 
assignments are shown in Fig. O, in units of the 
hadronic radius tq and in physical units on the left 
and right margins, respectively. 
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The masses of the two lowest-lying glueballs 
obtained on anisotropic lattices can now be 
compared to results using more conventional 
techniques. A meaningful comparison can be made 
by expressing results from different simulations in 
units of a common scale, for which I have again 
chosen the had r onic radius rg. Thus, the results 
from refs. |128, 131, 132| have been expressed or 
converted in units of ro, and whenever necessary 
the continuum extrapolation has been re-done. 
The resulting masses for the lowest scalar and 
tensor glueballs are listed in Table || together with 
the value of the aspect ratio f and the year in 
which the calculation was carried out. It is also 
indicated whether a continuum extrapolation has 
been performed. 

The table shows that the results obtained in 
different simulations are in broad agreement. If 
those estimates are to be used for phenomenological 
purposes, it should be kept in mind that they are 
valid in the quenched approximation. This implies 
that the conversion into physical units (e.g. by 
using rg — 0.5 fm) is ambiguous. Also, the influence 
of dynamical quark effects, and, perhaps most 
importantly, the issue of glueball-meson mixing has 
not been addressed for the results presented in 
Tabl e ^ . The latter has been studi ed, e .g. in 



also contains a detailed calculation of the spectrum 
of light quarkonia in quenched QCD. Overall, it 
is found that physical glueball masses are quite 
sensitive to the details of the proposed mixing 
pattern. 

Dynamical quark effects in th e glueball 
spectrum hav e b een studied in [137| and p8|. 



ref. 1 133 -135 1, and very recently in ref. |13f:], which 



Whereas ref. [137| reports no significant deviation 
of glueball masses computed for nf = 2 (within 
large errors), results by UKQCD Q indicate much 
lower estimates in unquenched simulations for both 
the scalar glueball and quarkonium state. However, 
a number of effects, e.g. the question of lattice 
artefacts have to be addressed in much more detail 
before these results can be confirmed. 

5.2. Heavy quarkonia and hybrids 

In addition to glueballs the spectrum of quarkonia 
and hybrids can also be calculated in lattice QCD. 
Such calculations may provide hints especially 
for the heavy quark sector, where very little 
experimental data exists for hybrid states. Thus, 
on the lattice one would like to compute the masses 
of states which are obtained from quark bilinears by 
inserting one or more gluons, such as bgb and cgc. 

Ideally a fully relativistic treatment of the heavy 
quarks would be desirable. However, for the 
currently accessible range of lattice spacings, for 
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Table 3. Comparison of the two lowest glueball masses in units of tq. 



Collab. 



ro TO2++ 



^ a — > year 



M+p 

GFll IM2 
Teper |28[ 
UKQCD 



4.21(11)(4) 5.85(2)(6) 3,5 V 1999 

4.33(10) 6.04(18) 1 Ij 1999 

4.35(11) 6.18(21) 1 V 1998 

4.05(16) 5.84(18) 1 x 1993 



which a ^ 0.05 fm (corresponding to ^ 4GeV) 
one expects large lattice artefacts in the charm 
sector, while the h quark with a mass above 4 GeV 
cannot at all be simulated directly. If anisotropic 
lattices with coarse spatial lattice spacings of — 
0.2 — 0.4 fm are employed, the situation is even 
worse. 

One way to address the problem of large cutoff 
effects in the heavy quark sector is to reso rt to 
an effective, non-relativistic treatment |13S], by 
introducing a cutoff A such that 



4.5 



(40) 



where mg is the mass of the heavy quark. In other 
words, relativistic states above A are excluded. 
Based on this approximation one can write down 
a discretised, effective, non-relativistic QCD action 
(NRQCD). On a formal level this approximation 
of QCD can be viewed as an expansion in the 
strong coupling constant and the 4-velocity v of 
the heavy quark. In this formulation the lattice 
spacing acts both as the UV regulator and the 
non-relativistic cutoff. This implies in turn that 
the continuum limit a — )■ cannot be taken, since 
otherwise the non-relativistic approximation breaks 
down. Therefore, in order for cutoff effects to be 
under control one relies on a "window" in a where 
NRQCD works well and lattice artefacts are small 
at the same time. 

Two groups have recently reported results for 
heavy quarkonia a nd hybri ds using NRQCD on 
anisotropic lattices [ 139 140{ . Both have used the 
O(a^) mean-field improved gluon action for ^ = 3, 5 
and an NRQCD action expanded in the 4-velocity 
to order rnqv^. At this order, spin interactions are 
neglected. Therefore one expects spin degeneracies 
for the S- and P-wave quarkonium states, as well 
as for each of the hybrid states Hi , H2 and H3 with 
quantum numbers J^^ = 1 (Hi), 1+"*" {H2) and 
0+"'' {H3), respectively. 

A convenient way to present the results is to 
quote the hybrid-5'-wave splittings S{Hi — S),i = 
1,2,3. In Fig. O those splittings, normalised to the 
IS — IP quarkonium splitting, are plotted versus 
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Figure 14. The hybrid splittings in the charm and 
bottom sectors from NRQCD on anisotropic lattices. 
Diamonds denote the results using relativistic heavy 
quarks. 



the spatial lattice spacing as. Focussing on the 
upper part of the figure which shows the spectrum 
in the bottomonium sector, one sees that anisotropy 
effects appear to be under control: the results by 
CP-PACS (l39| obtained for ^ = 3, 5 are consistent 
within errors. Furthermore, the data from both 
groups support the existence of a window for as ~ 
0.1 — 0.2 fm, where discretisation effects are small: 
in that range the variation of the results with 
is roughly as large as the statistical precision. It 
is obvious, though, that this is no longer the case 
for as > 0.2 fm. Similar observations apply in 
the charmonium sector shown in the lower part of 
Fig. |l^, where the scaling window appears to be 
somewhat larger. 

The results for the lowest bgb hybrid splitting 
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obtained on anisotropic lattices are [139, 140| 
1.542(8) GeV CP-PACS 



6{H, - S) 



1.49(2) (5) GeV KJM 



(41) 



For both charmonium and bottomonium hybrids 
the results can be compared to those in which the 
heavy quarks were treated relativistically, invol ving 
extr apolations in the heavy quark mass 141 , 
|l42|] . These data points are included as diamonds 
in Fig. |l^ and are consistent with those using 
anisotropic actions. This nicely illustrates how 
complementary formulations of heavy quarks can 
be used to control different systematic effects. 

6. Other topics 

The topics which are presented in this section could 
not be reviewed extensively, but nevertheless I shall 
briefly summarise the current status and refer the 
interested reader to the recent literature. 

6.1. Kaon weak matrix elements 

Lattice calculations for the _B-parameter Bk , which 

is relevant for K*^ — K mixing, as well as matrix 
elements for K — > tttt ande'/e have recentl y bee n 
reviewed by Kuramashi [ [143[ and Martinelli |144|. 

There are now many calculations of Bk using 
different discretisations of the Dirac operator. For 
staggered fermions the most recent result quoted in 
the naive dimensional reduction (NDR) scheme is 



rNDR 



(2 GeV) = 0.628 ±0.048. 



(42) 



Further improvements of this result could be 
achieved through the implementation of non- 
perturbative rcnormalisation. Simulations using 
Wilson fermions [145| are consistent with the result 
in eq. (42). Estimates for Bk using Domain Wall 
fermionsnave also been presented p^,|47[|. 

Compared to K'^ — K mixing the lattice 
predictions for matrix elements relevant for K — > 
TTTT and e'/e are a lot less accurate, and a number 
of systematic effects have to be b etter controlled. 
Details can be found in |143, 144 1. A surprising, 
negative result for Re(e'/e) has been reported 
recently in |4^. However, given that it has been 
obtained using the relatively new technology of 
Domain Wall fermions, future studies are required 
to clarify this issue. 



6.2. B decay matrix elements 

Decays of heavy-light mesons have long been 
studied in lattice gauge theories. The current status 



has been reviewed by Hashimoto ]146| |. Scmi- 
leptonic heavy-to-light and heavy-to-heavy decays 
have been analysed in the quenched approximation 
for several formulations of heavy quarks. Lattice 
estimates for heavy-light decay constants such as 
/b have stabilised in the quenched approximation, 
and their current values can be summarised as [146| 



fs = 170±20MeV, fsJIs = 1.15 ±0.04. (43) 

Present activities center on the quantification of sea 
quark effects. Preliminary estimates suggest that 
heavy-light decay constants increase by up to 20 % 
for nf = 2 flavours of dvnamica l qu arks at non- 
zero values of the lattice spacing |147|. However, a 
detailed scaling analysis is still required in order to 
separate sea quark effects from lattice artefacts, so 
that the large increase in /b for Uf = 2 does, in my 
view, not constitute a solid result at present. 

6.3. Structure functions 

This topic has recently received a lo t of attention 
and has been reviewed by Petronzio |145]. Current 
activities include the implementation of non- 
perturbative rcnormalisation of parton density 
operators, using either the RI or SF schemes 
described in section ^. The SF scheme has been 
employed to determine non-perturbatively the scale 
dependence of t he t wist-two, non-singlet parton 
density operator |14£]. Other projects in this area 



include calculations of highe r tw ist contributions to 
the pion structure function |150|. 



7. Summary 

The most significant theoretical development in 
Lattice Gauge Theory has surely been the progress 
made in formulating chiral symmetry at non-zero 
lattice spacing. It is now clear that chiral gauge 
theories can be put on the lattice in a consistent 
way, without breaking the gauge symmetry. This 
shows that a regularisation of the Standard Model 
exists beyond perturbation theory. 

Lattice simulations of QCD are becoming ever 
more refined, thanks to a number of technical 
developments, such as the implementation of 
the Symanzik improvement programme, non- 
perturbative rcnormalisation, the use of anisotropic 
lattices and bigger, more efficient simulations with 
dynamical quarks. 

The quenched approximation, which is still 
widely used, works surprisingly well. This is good 
news for many computationally more demanding 
applications, for which the quenched approximation 
will be useful for some time in the future. There 
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are now many attempts to quantify the effects of 
dynamical quarks, and in some cases they have been 
found to be significant. This has only been possible 

after all other systematic effects, including lattice 
artefacts, could be controlled at the level of a few 
percent. 

Although lattice QCD may not yet be the 
ultimate phenomenological tool, it is clear that 
enormous progress towards this goal has been, and 
will be made. 
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